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i2.  SPONSORING  military  activity 

Department  of  the  Army 

Office  of  the  Chief  of  Engineers 
Washington,  D.  C.  2031k 

In  tlMis  report  an  elasto-plastic  solution  is  given  for  determining  the  stresses, 
strains,  and  displacements  aroimd  a  circular  tunnel  in  rock.  The  relationship 
between  the  principal  stresses  in  the  plastic  zone  are  determined  by  the  Couiomb- 
Navler  failure  criteria  and  the  concept  of  "normality"  is  used  to  relate  the  plastic 
strains  to  the  "yield'  surface.  The  "normeility"  concept  results  in  a  dllatancy 
of  the  material  at  yielding  when  the  angle  of  shearing  resistance  of  the  material 
is  greater  than  zero.  Thus  the  displacements  calculated  by  the  method  given  in  this 
report  include  the  effect  of  dllatancy  In  the  yielded  zone  around  the  tunnel. 

'The  results  of  the  study  are  presented  in  graphical  form  which  enable  a  wide 
variety  of  problems  to  be  solved  by  hand  without  the  use  of  a  digital  computer. 

Example  problems  are  given  In  which  the  method  is  applied  to  various  situations  includ¬ 
ing  lined  tunnels  in  a  rock  mass  which  contains  a  destressed  or  loosened  zone  around 
the  tunnel  before  installation  of  the  liner.  The  method  of  analysis  presented  is 
especially  applicable  to  problems  where  loads  are  imposed  after  the  liner  is  installed 
as  in  the  case  of  the  design  of  protective  structures.  The  method  given  herein  is 
useful  for  design  because  the  calculation  of  deformations  and  strains,  rather  than 
stresses,  enable  the  designer  to  determine  if  a  particular  type  of  liner  can  tolerate 
the  deformations  imposed  at  the  rock-liner  Interface. 

Suggestions  are  also  given  for  the  selection  of  rock  mass  properties  to  be  lu^ed 
in  the  analysis,  i.e.  deformation  moduli  and  shear  strength  properties. 
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In  this  report  an  elasto-plastic  solution  is  given  for  determining  the 
stresses,  strains,  and  displacements  around  a  circular  tunnel  in  rock.  The 
relationship  between  the  principal  stresses  in  the  plastic  zone  are  determined 
by  the  Coulomb-Navi er  failure  criteria  and  the  concept  of  "normality"  is  used 
to  relate  the  plastic  strains  to  the  "yield"  surface.  The  "normality"  concept 
results  in  a  dilatancy  of  the  material  at  yielding  when  the  angle  of  shearing 
resistance  of  the  material  is  gi-eater  than  zero.  Thus  the  displacements  calcu¬ 
lated  by  th.'  method  given  in  this  report  include  the  effect  of  dilatancy  in  the 
yielded  zone  around  the  tunnel. 

The  results  of  the  study  are  presented  in  graphical  form  which  enable 
a  wide  variety  of  problems  to  be  solved  by  hand  without  the  use  of  a  digital 
computer.  E.xample  problems  are  given  in  which  the  method  is  applied  to  various 
situations  including  lined  tunnels  in  a  rock  mass  which  contains  a  destressed  or 
loosened  zone  around  the  tunnel  before  installation  of  the  liner.  The  method  of 
analysis  presented  is  especially  applicable  to  problems  where  loads  are  imposed 
after  the  liner  is  installed  as  in  the  case  of  the  design  of  protective 
structures.  The  method  given  herein  is  useful  for  design  because  the  calcxila- 
tion  of  deformations  and  strains,  rather  than  stresses,  enable  the  designer  to 
determine  if  a  particular  type  of  liner  can  tolerate  the  deformations  imposed 
at  the  rock-liner  interface. 

Suggestions  are  also  given  for  the  selection  of  rock  mass  properties  to 
be  used  in  the  analysis,  i.e.  deformation  moduli  and  shear  strength  properties. 
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NOTATION 


E  . 
seis 


=  Radius  of  a  circular  tunnel  or  tunnel  lining 
=  Constant  of  integration 
=  Tunnel  diameter 

=  Young's  modulus  of  an  elastic  medium  under  plane  strain  conditions 
=  Young's  modulus  of  a  loosened  or  destressed  zone  around  a  tunnel 
=  Young's  modulus  for  an  elastic  material  under  plane  stress  conditions 
=  Effective  Young's  modulus  of  deformation  for  a  rock  mass 
=  Young's  modulus  for  steel 

=  Dynamic  value  of  Young's  modulus  calculated  from  seismic  surveys 

=  Thickness  of  steel  lining 

1  sin  d 
~  1  -  sin 

=  Difference  between  the  circumferential  and  radial  stresses  in  the 
plastic  zone,  (Og  -  o^) 

=  Radial  pressure  on  the  surface  of  a  rock  tunnel 
=  Uniform  free-field  stress 
=  p  +  T 

=  Unconfined  compressive  strength 
=  Radius  to  the  elastic-plastic  boundary 
=  Radius  of  loosened  or  destressed  zone 

=  Radius  of  elastic-plastic  boundary  in  the  loosened  medium,  if  any, 
when  R  >  R^ 

=  Radial  distance  from  the  center  of  a  tujinel 
=  Stress  difference  in  the  elastic  region  (og  -  o^) 

=  Joint  spacing 


N  -  1 


Radial  displacement 
Tangential  strain 
Tangential  plastic  strain 
Radial  strain 
Radial  plastic  strain 

Poisson's  ratio  of  medium  in  plane  strain 
Poisson's  ratio  of  loosen*.!  or  destressed  raeiium 
Poisson's  ratio  for  plane  stress 
Angle  of  shearing  resistance  of  the  material 
Radial  stress 

Radial  pressure  axerced  by  steel  liner  against  the  adjacent  rock 
medi urn 

Circumferential,  stress  or  tangential  stresses 
Radial  stress  at  radius  R 

Unconfined  compressive  strength  of  the  rock  mass  surroimding  a  tunnel 

Yield  stress  of  steel 

Shear  stress  in  polar  co-ordinates 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 


British  units  of  measurement  used  in  this  report  cun  be  converted  to  metric 


units  as  follovs ; 

Mult inly 

BX. 

To  Obtain 

inches 

2.5** 

centimeters 

feet 

0.30U8 

meters 

cubic  Inches 

16.3871 

cubic  centimeters 

pounds 

0.*(5359237 

kilograms 

pounds  per  squeu-e  inch 

0.070307 

kilograms  per  square 
centimeter 

pounds  per  cubic  foot 

16.0185 

kilograms  i>er  cubic 
meter 

inch-pounds 
inches  per  second 


0.011521 

2.5^ 


meter-kilograms 
centimeters  per  second 


Chapter  1 


Introduction 


General 


Cne  of  the  most  significant  engineering  properties  of  rock  masses  is  shear 
strength.  The  stability  of  both  lined  and  unlined  openings  under  various  static 
and  dynamic  loadings  depends  on  the  shear  strength  of  the  surrounding  mass. 

Since  the  shear  strength  of  rock  masses  increases  with  confining  pressure,  it  is 
necessary  to  employ  an  analysis  which  includes  this  property  in  order  to  properly 
assess  the  effects  of  various  types  of  tunnel  linings. 

One  of  the  first  attempts  to  take  into  account  the  Coulomb-Navier  shear 
strength  characteristics  of  earth  materials  in  a  yielded  zone  around  a  deep 
tunnel  was  made  by  Terzaghi  (1919)*  Although  he  did  not  obtain  a  general  solu¬ 
tion  to  this  problem  he  was  impressed  by  his  field  observation  that  deep  bore¬ 
holes  in  frictional  materials  remained  stable  at  depths  where  the  material 
adjacent  to  the  opening  should  have  failed.  His  curiosity  about  this  problem 
led  him  to  ask  H.  M.  Westevgaard,  a  colleague  of  his  at  Harvard,  two  questions. 

(1)  "What  distributions  of  stress  are  possible  in  the  soil  around  an 
unlined  drill  hole  for  a  deep  well?" 

(2)  "What  distributions  of  stress  make  it  possible  for  the  hole  not  to 
collapse  but  remain  stable  for  some  time  either  with  no  lining  or 
with  a  thin  "stove-pipe"  lining  of  small  structural  strength?" 

Westergaard  (19^0)  published  the  answers  to  these  questions  after  obtaining  a 
solution  for  tne  stresses  around  a  borehole  where  the  material  in  the  plastic 
zone  around  the  hole  was  assumed  to  fail  according  to  the  Coulomb-Navier 
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failure  criteria.  The  results  showed  that  a  small  radial  confining  pressure  at 
the  surface  of  the  borehole  enabled  the  radial  stresses  to  increase  rapidly  with 
depth  behind  the  surface  such  that  the  radial  pressures  a  few  inches  behind  the 
borehole  surface  furnished  sufficient  confinement  to  support  the  high  circum¬ 
ferential  stresses  around  the  opening.  The  results  of  Westergaard' s  solution  were 
expa.nded  and  interpreted  by  Terzaghi,  19^3.  Essentially  the  same  solution  has  been 
published  by  Jaeger  (1956),  Jaeger  and  Cook  (I969)  and  Sirieys  (196I4).  In  this 
solution,  the  stresses  are  obtained  around  a  long  circular  tunnel  in  a  medium 
where  the  free-field  principal  stresses  are  equal  (i.e.  =  Ojj  =  p^.  Fig.  l). 

For  this  case,  the  shear  stresses  around  the  opening,  are  equal  to  zero  and 

the  differential  equation  of  equilibrium  shown  in  Fig.  1  can  be  expressed  in  terms 
of  the  radial  stress  0  and  circumferential  stress  0^.  The  -olution  was  obtained 
by  integrating  the  differential  equation  of  equilibrium  in  the  plastic  zone  (Fig. 

1)  with  the  constraint  that  the  circumferential  and  radial  stresses  are  related 
by  the  Coulomb-Navier  failure  criteria  and  requiring  that  the  radial  stresses  are 
continuous  at  the  boundary  between  the  elastic  and  plastic  zones  (Fig.  l) . 

An  illustration  of  the  stress  distributions  obtained  from  the  Westergaard 
solution  is  given  in  Pig.  2  for  a  tunnel  in  sand  (C  =  0,  (J>  =  30°).  The  distribu¬ 
tions  of  circumferential  stress  and  radial  stress  o  are  shown  in  dimensionless 

form  in  terms  of  the  free-field  stress  ana  three  differenL  stress  distributions 
are  shown  for  cases  where  the  confining  press'ore  on  the  inside  of  the  tunnel  is 
equal  to  1/10,  1/20,  and  1/1*0  of  the  free-field  stress  p^.  For  each  stress  dis¬ 
tribution  shown  the  circumferential  stress  is  maximum  at  the  boundary-  of  the 
elastic  and  plastic  zones  and  the  distance  to  the  boundary  between  the  elastic 
and  plastic  zones  increases  as  the  ratio  of  the  free-field  stress,  p^,  to  the 
internal  pressure,  p^ ,  increase.s.  It  should  also  be  noted  that  the  circumferential 
stress  increases  very  rapidly  with  depth  behind  the  opening  in  the  plastic  zone 
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and  illustrates  the  atility  of  a  frictional  material,  to  carry  high  circumferential 
stresses  at  depth  due  to  providing  a  nominsl.  confining  pressure,  p^ ,  on  the  inside 
surface  of  the  tunnel. 

The  solution  discussed  above  however  yields  only  the  elastic-plastic  stress 
distribution  around  the  opening.  A  knowledge  of  the  stress  distribution  is  only  of 
academic  interest  for  the  design  of  a  tunnel  lining  unless  the  deformations  at 
the  wall  of  the  tunnel  axe  compatible  with  the  deformations  a  lining  can  resist 
before  failing.  The  radial  displacements  at  a  tunnel  wall  are  due  to  both  the 
inelastic  strains  in  the  plastic  zone  adjacent  to  the  tunnel  eind  the  elastic 
strains  in  the  rock  outside  the  plastic  zone.  Although  the  strains  in  the  elastic 
region  are  easily  calculated,  the  calculation  of  strains  in  the  plastic  zone 
involves  some  assumption  regarding  the  relations  between  the  plastic  strains  and 
volume  changes  of  the  material  at  failure.  Newmark  (1969)  obtained  a  closed 
solution  for  the  problem  described  above  and  shown  in  Fig.  1  which  considers  the 
displacements  and  strains  around  the  opening  as  well  as  the  stress  distributions. 

To  obtain  this  solution  Newmark  combined  the  differential  equation  of  equilibrium 
and  the  compatibility  equation  relating  radial  displacements  to  radial  strains  and 
circumferential  strains  to  yield  the  relation  given  as  equation  (U)  in  this  report. 
Since  the  differential  equation  of  equilibrium  and  the  compatibility  relation 
involve  only  equilibrium  and  geometry  respectively  it  follows  that  the  relation 
given  by  equation  (U)  must  be  satisfied  in  both  the  elastic  and  plastic  regions. 
Newmaik  (I969)  obtained  a  solution  for  equation  (U)  by  assuming  that  after  yield 
of  the  material  that  the  plastic  components  of  the  radial  and  circumferentia.1 
strains  were  such  that  the  resulting  volume  change  was  zero  or  the  yielded  material 
behaved  as  an  incompressible  solid.  The  results  of  the  solution  have  been  pre- 
paxed  by  Newmark  (1969)  in  the  form  of  charts  from  which  the  displacements  at  the 
tunnel  wall  can  be  determined  as  a  function  of  the  elastic  properties  of  the 
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medixun  (v,  E),  the  shear  strength  properties  of  the  medim  as  determined  by  the 
unconflned  compressive  strength  and  the  angle  of  internal  friction,  the  free- 
field  stress,  and  the  capacity  of  the  tunnel  lining.  Since  real  rock  masses 
increase  in  volume  at  failure,  a  phenomenon  called  "dilatancy",  it  is  probable 
that  the  Newmark  solution  underestimates  the  inward  radial  displacements  at  the 
tunnel  wall  because  of  the  assumption  that  the  plastic  components  of  strain  yield 
no  vol’ome  change. 

Scope 


In  this  report  an  elasto-plastic  solution  is  given  in  Chapter  2  which 
accounts  for  the  dilatant  properties  of  a  rock  mass  which  obeys  the  Coulomb- 
Navier  failure  criterion.  Example  calculations  are  given  to  illustrate  the 
variety  of  problems  which  may  be  solved  by  this  method.  Suggestions  for  the 
selection  of  rock  mass  properties  and  conclxisions  are  given  in  Chapter  3. 
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Chapter  2 
Analysts 

Tlie  basic  configuration  considered  in  this  study  is  shown  in  Fig.  1.  This  illus¬ 
tration  represents  a  section  throiigh  an  infinitely  long  tiinnel  where  the  strain 
in  the  direction  of  the  tunnel  axis  is  zero.  The  cylindrical  coordinate  system 
with  the  origin  at  the  center  of  the  circular:  opening  is  used  in  this  analysis. 


Fundamental  Relations  of  Equilibrium  and  Compatibility 

The  relations  given  herein  are  derived  for  plane  strain  conditions  for  a 
uniform  stress  field. 

The  basic  differential  equation  of  equilibrium  for  a  typical  element 
shown  in  Fig.  1(a)  is 


0 


(1) 


where  and  Og  represent  the  radial  and  tangential  stresses  respectively,  at  a 
radial  distance  r. 

The  radial  and  tangential  strains  for  the  assumed  conditions  of  pleuie 
strain  can  be  stated  in  terms  of  the  radial  displacement,  u,  as  follows: 


z 

r 


9u 

3r 


= 


u 

r 


(2) 


Upon  eliminating  the  displacement  u  from  Eqs.  (2),  one  obtains  the 


following  equation  of  compatibility: 


Eqs.  (l)  euid  (3)  can  now  be  combir.ed  into  a  single  equation  by  multiplying  Eq.  (l) 
by  the  quantity  -  (1  -  y)  r  and  adding  it  to  Eq.  (3)  multiplied  by  the  quantity  E. 
The  result  is : 

It  is  to  be  noted  here  that  Eqs.  (l),  (3)  and  (It)  are  valid  for  both  elastic  and 
inelastic  conditions.  However,  the  validity  of  Eqs.  (?),  (3)  and  (it)  is  limited 
to  the  case  of  small  displacements. 

In  order  to  solve  Eq.  (it),  it  is  necessary  to  relate  the  term  in  brackets 
at  the  extreme  right  of  Eq.  (it)  uo  the  other  terms  in  that  equation  through  the 
use  of  stress-strain  relationships  and  to  establish  the  appropriate  boundary 
conditions . 

Relations  for  the  Elastic  Zone 

The  stress-strain  relationships  for  the  elastic  behavior  of  a  material 
under  plane  strain  conditions  can  be  written  as  follows; 


Ee^  =  Cq  - 

(5) 

ECr  =  0^  -  pog 

(6) 

where  E  and  p  are  the  Young's  modulus  and  Poisson's  ratio  of  the  material  under 
plane  strain  conditions.  I'he  values  of  E  and  p  for  plane  strain  are  related  to 
the  corresponding  values  for  plane  stress  by  the  following  equations : 

E  =  E  /(I  -  p  ^)  (7) 

o  o 

p  =  p  /(I  -  p  )  (8) 

0  o 
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where  E  is  the  Young's  modulus  for  plane  stress  and  u  is  the  Poisson's  ratio 
o  o 

for  plane  stress.  If  the  quantity  S,  as  suggested  by  Newmark  (1969) »  is  defined 

S  =  EEq  -  (1  -  p)  (9) 

then  from  Eq.  (5)  it  follows  that 

S  =  (o^  -  o^)  (10) 

Thus  in  the  elastic  region  the  quantity  3  is  equal  to  the  stress  difference 
between  the  tangential  stress  o.  and  the  radial  stress  o  .  Eq.  (6)  gives 

Ee^  -  (1  -  p)  Og  =  =  -  s  (ll) 

Using  Eqs.  (9)  and  (ll),  Eq.  (4)  can  be  simplified  as; 

r  ^  +  2S  =  0  (1£) 

+  2  (o.  -  a  )  =  0  (ISi 

9  r 

2  (o^  -  Oq)  =  0  (14) 


or 


r  -r-  (o  -  o  ) 
or  o  r 


Fi-om  Eq ,  ( 1 ) 


3o 

2r  + 

or 


Adding  Eqs,  (13)  and  (l4). 


Eq.  (15)  indicates  that  in  the  elastic  region,  the  sura  (a.  +  a  )  is  a  constant. 

D  r 

At  very  large  values  of  r,  each  of  the  stresses  and  is  equal  to  the  free- 
field  uniform  stress,  p^.  Therefore  in  the  elastic  region. 
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(16) 


Relations  for  the  Plastic  Zone 

The  radial  and  circuniferential  (or  tangential)  strains  in  the  plastic 
region  can  be  vritten  as  the  sum  of  the  respective  elastic  and  plastic  components. 
Thus  in  the  plastic  zone, 

=  Zq  (elastic)  +  (plastic) 

(21) 

and  e  =  c  (elastic)  +  e  .  (plastic) 

T  T  1 
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The  concept  of  perfect  plasticity  requires  that  the  associated  strain  rate  vector 
must  be  normal  to  the  yield  surface,  Drucker  and  Prager  (1953).  Accordingly,  if 
f  represents  the  yield  function  which  is  valid  for  the  plastic  zone,  the  plastic 
strain  rate  components  in  the  radial  and  tangential  directions  can  be  related  by 
the  equation 


Cq  (plastic)  af/SOg 

I  (plastic)  ~  af '/  Do 
r  r 


(22) 


The  yield  function  used  in  this  analysis  is  based  on  Coulomb-Navier  yield  criterion 
and  is  given  by 


f  =  C,  -  ”  0  (23) 

1  +  s  n  6 

where  ^  =  tinconfined  compressive  strength  of  the  material 

and  4i  =  angle  of  shearing  resistance  of  the  material.  Thus  Eq.  (22)  gives 


Cq  (plastic) 
(plastic) 


constant . 


Since  the  ratio  of  the  plastic  strain  rate  components  is  a  constant  (for  the  yield 
function  used  in  this  analysis),  the  ratio  of  the  plastic  strain  components  will 
also  be  equal  to  the  same  constant.  Thus  the  ratio  between  the  plastic  strain 
components  may  be  written  as 


£g  (plastic) 

1 

(plastic) 

or 

Eq  (plastic) 

=  "ep  =  ~  JT  '  S  (plastic) 

or 

(plastic) 

=  -  N  •  F 

4>  '6p 

(2h) 
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It  should  be  noted  that  the  relation  given  in  Eq.  (2k)  produces  an  Increase  in 
volume  at  feiilure  and  the  percentage  of  volume  change  increases  as  N.  increases. 
Eq.  accounts  for  the  primary  difference  between  the  solution  presented  herein 

and  the  solution  given  by  Newmark  (1969)  where  the  condition  of  incompressibility 
assiimed  by  Newmark  imposed  the  condition  that  =  0.  By  utilizing 

Eq.  {2k)  Eqs.  (21)  can  now  be  rewritten  as: 


^0 

=  1  (<^0  -  w<^r)  + 

^0p 

(25) 

c 

r 

=  1  (0^  -  ncg)  - 

%  ^9p 

(26) 

=  Og  -  tja^  +  E  • 

Sp 

(27) 

Ee 

r 

=  "r  -  ‘“’e  -  % 

‘  *  ^ep 

(28) 

Substituting  Eqs.  (27)  and  (28)  into  Eq.  (U) 

r  If  .  2S  .  E  Cep  -  1)  =  0  (29) 

where  S  =  EEq  -  (l  -  p) 

The  quantity  S  is  equal  to  the  difference  between  the  tangential  6uid  radial  stresses 

in  the  elastic  range  but  in  the  inelastic  (plastic)  range  it  has  only  a  formal 

meaning.  For  the  inelastic  conditions,  the  difference  between  the  tangential 

stress  and  the  radial  stress  o  is  defined  as: 

0  r 

P  =  Og  -  o^  (30) 

By  Eq.  (27).  Eeg^  =  Ee^  "  *^6  ^'^r  “  ^  ~  °r^  ~  S  -  P  (3l) 
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'i'hereforo  Eq.  (29)  becomes 


r  ■—  +  2S  +  (S  -  F)  (tl  -  1)  =  0 

dr  4 


or 


>■  S  *  S  (N  t  1) 


P  (B,  -  1) 


(32) 


It  may  be  noted  that  Eq.  (32)  is  valid  in  tae  plastic  region  only. 

In  the  plastic  region  the  Coulomb-I/avier  yield  criterion  is  assumed  to 
be  valid.  Accordingly 


o„  =  0  •  N.  +  0 

6  r  if  u 


(33) 


P  =  o.  -  a  =  (N.  -  1)  o  +0 

e  r  (J)  r  u 


(3U) 


Differentiating  Eq.  (3^0  with  respect  to  r 


aP 

=  (N  -  1) 

3r  '  <t  3r 


(35) 


Substituting  Eq.  (35)  into  Eq.  (1) 


r  -  (N^  -  1)  P  -  0 

dr  <fc 


(36) 


The  solution  of  Eq.  (36)  yields 


P  =  P  - 
a  a 


N^-1 


whore  P  =  P 
a 


r=a 


(37) 
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Let  the  radial  pressure  inside  the  tunnel  be  and  the  radius  of  the  circular 
tunnel  be  a.  Then 


Therefore 


P  =  P 
a 


p  s  p 


r=a 


=  ^  ‘  Pi 


■  iNa-i  r  -I  r  ^ 

— I  a  O  +  (N  -  l)  p  [— 

aj  [_  u  (fi  *^i  [aj 


Na-1 


(38) 


(39) 


At  any  radius  r  within  the  plastic  zone, 


°r  “  (P  -  au)/(H  -  1)  = 


•^u  1  frl^-^  % 


i  N ,  -  1  al  Ua  -  1 


=  (Pi  +  T) 


Na-1 


-  T 


where  T  =  o 


■  N.  •  o  +  o 
«  9  r  u 


/(ir  -  1) 

u'  (i 


and 


(UO) 

(kl) 

(33) 


Substitution  of  Eq.  (39)  into  Eq.  (32)  gives 


^  H  2  -  1) 

3r  9 


[°u  ^  Pi  [i 


Na-1 


N*-l 


where 


A  = 


=  A  r 


(N,  -  l)'^  (p.  +  T) 
<P  1 


Solution  of  Eq.  {U2)  is  of  the  form 


N,+l 

S  r 


2N, 


2Na 

r  ^  +  C 


(42) 


(43) 


(44) 
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vhoi'e  C  is  the  coiistojit  of  integration.  Tlie  constant  C  has  to  be  evaluated  from 
the  boundary  conditions  specified  for  any  given  case. 

Seven  specific  coses  have  been  studied  in  this  investigation,  each  with 
a  different  set  of  boundary  conditions.  Tliese  cases  are  presented  and  analyzed 
in  the  remainder  of  this  chapter. 


N.  -  1 


Pu  =  o  +  (N  -  1)  0  =  0  +  °  •  (2p  -  0  ) 

h  u  4  R  u  N,  +  1  ■^o  u 


=  ir-TT  'Fo  <“♦  - 


Equating  Eqs.  (45)  and  (4?) 


,N.-1 

r)  “P 


2  Pq 

♦  1  ■  p,  *  T 
f  1 


Po  T  N^-1 


P.  ^T| 


(N  -  1) 

At  r  =  R.  S,  =  =  2  (p^  T) 


Now  the  value  of  the  constant  in  Eq.  (44)  can  be  estimated. 


According  to  Eq.  (44),  at  r  =  R 


C  =  .  R 


A  ■  R 

2N, 


Substitution  of  Eq.  (50)  into  Eq.  (44) 


A  r 

2Na 


2N^  2N.1  l> 

r  ^-R  ^  +S„*R 

K 


Substituting  for  A  and  S„  in  Eq.  (51),  and  simplifyint^ 

K 


_ S__  r 

P,-  T  [a. 


-2- -  L  .(N  _  1)  +  £  (M  +  1)  (52) 

2N^  [aj  ^  (aj 


The  value  of  R/a  to  be  used  in  Eq.  (52)  is  given  by  Eq.  (U8).  Eq.  (52)  gives  the 
value  of  S  at  any  radius  r  in  the  medium. 

When  the  values  of  S  €ir.d  at  any  radius  are  known,  one  can  calculate  the  cir¬ 
cumferential  strain  e-  at  reuiiua  r  using  the  following  relation: 


=  I  js  t  (1  -  ,) 


The  tangential  plastic  strain  component  at  any  radius  r  in  the  plastic  zone  cem 
be  estimated  from  the  relationship: 


c  =  i  (S  -  P) 
P  E 


where  S  and  P  are  given  by  Eqs .  (9)  and  (39)  respectively. 

Eqs.  (52)  and  (39)  can  be  combined  into  a  single  equation  to  yield  the  S  -  P 
relationship  applicable  in  the  plastic  region: 


N,  -  1 


N,  +  1 


R>  ^ 


The  corresponding  relationship  applicable  to  the  elastic  region  can  be  readily 
derived  from  Eqs.  (IT)  and  (37): 


1=  (f) 


N^+1 
i  V 


Eq.  (^*7)  yields 


I>0  * 


M.  ♦  1 

2  (n  -  1] 

9 
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Let  U3  now  define 

Q  = 

p  +  T 

0 

(58) 

(N.  +  1) 

and 

K  = 

-  >  -  —  - 

2  (N.  -  1) 

(59) 

'P 

Th\;is  Eq.  (57)  becomes 

Q  = 

According  to  Eq.  (37) 


It  is  to  te  noted  that  the  quantity  Q  is  independent  of  the  radial  stress 
whether  the  stress  sit\iation  is  plastic  or  elastic. 


Graphical  Solution 


Althoiigh  the  various  quantities  of  interest  for  design  can  be  calculated 
from  the  relationships  given  above,  charts  summarizing  the  relationship  between 
P,  S,  and  Q  can  be  made  to  facilitate  calculations,  Newmark  (1969).  Examples  of 
such  cherts  are  given  in  Figs.  3  through  7  for  values  of  ranging  from  2  to  6. 
These  charts  are  based  on  Eqs .  (55) >  (56)  and  (60). 
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The  abscissa  in  the  chart  is  S  and  the  ordinate  is  P,  both  being  plotted 
to  a  logarithmic  scale.  Lines  sloping  up  to  the  right  represent  constant  values 
of  R/r  and  the  lines  sloping  up  to  the  left  give  constant  values  of  Q.  The 
heavy  line  for  R/r  =  1.00  represents  the  limit  of  elastic  behavior.  Below  this 
line  the  behavior  is  plastic  and  above  this  line  the  behavior  is  elastic.  If  the 
circumferential  strain  and  the  radial  stress  are  known  at  any  point  with  radius  r, 
then  both  S  and  P  at  that  point  are  known  and  their  intersection  can  be  determined 
This  intersection,  if  below  the  heavy  line  for  R/r  =  1.00,  gives  immediately  the 
value  of  Q  from  which  may  be  determined.  The  value  of  R/r  is  also  obtained, 
which  can  be  used  directly  to  determine  the  radius  of  the  elasto-plastic  boundary. 
If  the  point  of  intersection  is  in  the  elastic  region,  the  original  veQue  of  P  is 
not  valid  b”t  can  be  determined  directly  from  Eq.  (l6a)  since  o^.  and  S  are 
known  to  begin  with. 

Similarly  if  the  free-field  stress  and  the  radial  stress  at  any  radius 
r  are  known,  the  quantities  Q  and  P  can  be  calculated  and  these  values  can  then 
be  used  to  define  the  initial  point  in  the  chart.  If  the  plotted  point  is  in 
the  plastic  region,  it  gives  immediately  the  values  of  R/r  and  S  from  which  the 
radiiis  of  the  elasto-plastic  boundary  and  the  circumferential  strain  at  radius 
r  can  be  determined.  If  the  plotted  point  is  in  the  elastic  region,  the  plotted 
value  of  P  is  not  valid  but  the  value  of  S  is  always  valid.  This  value  of  S  can 
then  be  used  to  determine  the  circumferential  strain  at  the  given  radius  r. 

If  we  now  proceed  from  a  given  radius  r  to  some  new  value  of  r,  say  r  ==  b , 

with  R  known,  then  R/b  can  be  calculated  and  this  value  read  os  the  new  value  of 
R/r.  One  proceeds  to  this  new  value  along  a  line  of  constant  Q  since  does  not 
change  in  the  same  material.  This  new  point  gives  directly  the  value  of  S  as  the 
abscissa  and  the  value  of  P  as  the  ordinate  for  the  new  radius.  If  the  inter¬ 
section  is  in  the  elastic  range  of  the  chart,  the  value  of  P  is  not  valid  but  the 
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value  of  S  is  alv&ys  vEilid.  If  the  intersection  is  in  the  plastic  region,  the 


new  P  is  valid  and  o  can  be  determined  directly  because  a  is  known.  With  o 

r  u  r 

determined,  e.  can  be  determined  directly  from  S.  If  the  new  P  is  not  valid,  or 
the  situation  is  elastic,  can  be  determined  from  Eq.  (l6a)  since  both  S  and 
are  known  for  the  new  radi\is .  Eq.  (53)  can  then  be  used  to  determine  the 
circumferential  strain  e„. 

Regardless  of  whether  the  starting  point  is  in  the  elastic  region  or  in 
the  plastic  region,  a  new  point  at  another  radius  can  be  determined  by  going  eilong 
a  line  of  constant  Q,  taking  into  account  of  the  fact  that  the  change  in  radius 
is  one  which  corresponds  to  a  constant  value  of  R,  which  can  be  determined  from 
the  first  point  located.  It  is  thus  apparent  that  one  can  start  at  any  point  and 
go  to  any  other  point  in  the  same  medium  directly  but  one  must  start  over  again  in 
each  new  medium,  with  some  valid  starting  point  being  known. 

The  preparation  of  a  chart  for  a  given  value  of  involves  a  large  number 
of  repetitive  calculations .  Thiis ,  a  small  computer  program  was  written  and  used 
to  generate  the  data  necessary  for  construction  of  the  cheu-ts  shown  in  Figs.  3 
throuigdi  7  which  were  prepared  for  values  of  =  2,  3,  5,  and  6.  If  the  value 

of  N  in  a  given  design  problem  is  not  an  integer,  the  calculations  may  be  done 
by  interpolation  using  the  appropriate  charts.  The  range  of  input  parameters 
covered  by  the  charts  shown  in  Figs.  2  through  6  is  large  enougli  to  include  nearly 
all  problems  involving  deep  protective  structui’es  in  rock. 
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*  =  37° 

Lki) 

(52) 

(1*8) 
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iliUDM^jU 


^e  ”  i"  Pi^ 


625.000 

6  X  10^ 


=  l.ou  X  10~^  or  10. U  % 
Graphical  Solution  (Use  Chart  for  Na  =  4) 


p  +  T 
^o 


r=a 


Q  =  17,067  psi 

+  3  •  =  2000  psi 


For  Q  =  17,067  psi  and  P  =  2000  psi,  read  from  the  chart  for  N,  = 

9 


^  =  2.17 


eind 


s  =  625,000 


e 


0 


r=a 


I  (S  +  (1-ii)  p.) 


=  10. %, 

I 

I 

Case  2.  Input  Parameters  are  the  Same  as  for  Case  1  except  that  bhe  value  of 

i 

0^  in  the  plastic  zone  varies  linearly  from  zero  at  the  tunnel  vail 
(r=a)  to  the  full  value  at  the  elasto-plastic  boundary  (r=R) 


The  assumption  of  a  linear  variation  of  in  the  plastic  region  intro¬ 
duces  the  following  changes  in  the  relationships  derived  for  the  plastic  region. 

The  Coulomb-Navier  yield  criterion  for  the  plastic  region  becomes  modified 
as 


a 


0 


a 

r 


+  a 

u 


(33a) 
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where  R  =  radius  of  elasto-plastic  boundary 

The  stress  difference,  P,  at  radius  r  in  the  plastic  region  is  now  given  by 


Substituting  Eq.  (35a)  into  Eq.  (l) 


The  solution  of  Eq.  (36a)  yields 


^  [aj  (r  -  a)(N.  -  2) 

'  9 


1^1 

a  "a 


(37a) 


where  P  =  P 

a 


=  p,  (»,  -  1) 


r=a 


-  “r  ■  »r  <"♦  -  ♦ 


u 


R  -  a 


(r  -  a) 


Differentiating  Eq.  (3^a)  \dLth  respect  to  r. 


3o 


u 


R  -  a 


C3l;a) 


C35a) 


1^  -  (N  -  1)  -  = 

3r  6  r  R  -  a 


(36a) 


Therefore 


W.-l 


P  =  Pi  (f)  ^  rR-'-VriJr-  2)’ 


At  any  radius  r  within  the  plastic  zone,  (a  ^  r  ^  R) 


N.-l 


^  -  ”u  •  F^j  / 


(N,  1) 


(39a) 


(liOa) 
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a.  =  a  •  N,  +  0  ^  — — 

e  r  (p  u  R  -  a 


(33a) 


Substitution  of  Eq,.  (39a)  into  Eq.  (32)  yields 


-  §  ^  S  (N^  .  1) 


N.-l 

A  r  ^  +  B  r 


(h2a) 


(N,  -  1)‘ 


where 


T  •  a 


euid  B  = 


N.-l  ^i  (R  -  a)tN,  -  2 


[R  -  a)  ("2  -  N  ,  : 


(43a) 


Solution  of  Eq.  (42a)  is  of  the  form 


2N.  N ,  +  2  ^ 

91  9 


(44a) 


where  C  is  the  constaint  of  integration. 
At  r  =  R,  S  =  Sp 


Therefore 


^R  • 


-L.  V  <t>  ^  B  •  R 
2N  ^  N.  +  2  ^ 

9)  <p 


Eliminating  C  from  the  above  equation,  Eq.  (44a)  can  be  rewritten  as 


S  •  r 


2N.  2N, 


-i-  r  -  R  B  r  9 


N,+2  N  +2 


+  Sp  .  R 


(51a) 
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At  the  elasto-plastic  boundary  (r=R) 


■■r  =  r-TT  (Po  * 


According  to  Eq.  (39a) 


P  =  p  (N  -  1)  7 

R  di  a 


Equating  Eqs.  (U7)  and  (lt5a) 


o  'a 
u 

(r  -  a)CN.  - 


(i<5a) 


Pi  fn  '  - 

la  4i 


2(p-  ♦  T) 

-n^rry 


k  I 


There  is  no  closed  solution  of  Eq.  (l»8a)  for  R/a.  Eq.  (l»8a)  has  to  be  solved  for 
R/a  only  by  successive  approximation. 


At  r  =  R, 


2  (P„  »  T)(»^  -  1) 

*  1) 

9 


Substituting  this  value  of  into  Eq.  (51a)  and  simplifying 


.  2N,  2N  _  N  +2  N.+2 

JL  T.  R  'J’  +  B  ^9  rj  9 

2N  ^  ~  ^  N  +  2  ■  ^ 

▼  I  J  V  i 


2  (p  +  T)(N.  -  1)  N>1 

+  - 2 - i -  .  R  9 

N.  +  1 
9 


(52a) 


where  A,  B,  and  R  are  given  by  Eqs.  (l»3a)  and  CU8a).  Eq.  (52a)  gives  the  value 
of  S  at  any  radivis  r  in  the  medium.  The  circumferential  strain  can  then  be 
calcialated  using  Eq.  (53). 


In  the  elastic  region  (r  ^  R)  the  stresses  and  strains  are  given  by 
Eqs.  (19).  (20),  (5),  (6)  and  (1*6). 

Example  2 

Data: 


Unlined  Tunnel: 


a 

P, 


u 


N, 


8* 

2000  psi 
6x10^  psi 
1/3 

1200  psi 
k 
0 


Required:  e 


=  ? 


r=a 


Solution: 


T  -  =  boo  psi 


(bl) 


boo 


5.1 

a 


R 

a 


boo 


5.  1 

a 


2  X  2b00 


(b8a) 


Solving  7  =  1.75 


or 


R  =  1.75  X  8  =  lb' 


Note :  If  0^  is  constant  over  the  plastic  region 


R 

a 


2  X  2b00l^^^ 

5  X  boo 


(b8) 


r  =  1.3b 


2b 


A  -  _2.  ^00 

^  J  0.75  3C  2 

EL 


1200  »  3 
a  •  0.75  •  2 


-  2U00 

a 


(U3a) 


2  X  21*00  X  3  _5 

+  _  R 


(52a) 


Sutatituting  —  =  1.75 
a 


=  300  [ 

1  -  1.75^ 

-  1*00 

1  -  1.75^ 

L 

r«a 

L 

»  J 

+  28800  X  1.75^ 

=  32200  psi 

-  i  [s  +  (1  -  p)  p^J 


i.l*  X  10~^ 


Note:  If  0^  ia  constant  over  the  plastic  zone 


=  1*00  X  -1 

3  +  5  X  1.31*®] 

=  8300  psi 

(52) 

r=a 

0  ^ 

=  i  (S  +  (1  -  p)  pj  =  1.1*  X  10~^ 
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Case  3.  Input  peirameters  are  the  same  aa  for  Case  1  except  that  cohesion  la 


asBumed  to  be  lost  in  the  entire  plastic  region,  i.e. , 

Ou  =  0 


The  assvanption  of  a  complete  loss  of  cohesion  in  the  plastic  region  results 
in  the  following  modified  equations: 

Coulomb'Navier  failure  criteria... 


o.  =  a  •  N. 
e  r  d 


Stress  difference  P  =■  a  (N.  -  l) 

r  4> 


(33b) 

(3Ub) 


P  a  P  . 
A 


rr)  ^ 


(37b) 


where 


P  a  P 
a 


Pi  ("♦  - 


r=a 


At  any  radius  r  within  the  plastic  region 


=  P/(N^  -  1) 

r  o 


°e  "  '^r  * 


(UOb) 

(33b^ 


At 


r  =  R,  Pp  =  N ,  +  1  Po 


(47b) 


=  Pi  f 


9 


:45b) 


Thus 


R 

a 


N,  +  1  p. 
d> 


9 


(48b) 
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Eq.  (52)  is  now  simplified  as 


N.+l 


2N, 


2N. 


2N, 


%  -  1)  + 


(52b) 


Eq.  (52b)  can  now  be  used  to  determine  the  tangential  strain  at  radius  r. 

In  the  elastic  region,  r  ^  R,  the  stresses  and  strains  eu*e  determined  using 

Eqs.  (19),  (20),  (5)  and  (6)  noting  that  o  =  2  p^/(N  +  l). 

n  09 


Case  I*.  Input  parameters  are  the  same  as  Case  1  except  that  the  quantities  E 
and  p  have  c'ifferent  values  in  the  elastic  and  plastic  regions. 


Let 

E,  E'  =  Yo^ing's  modulus  in  the  elastic  emd  plastic  regions 
jj,  p'  =  Poisson's  ratio  in  the  elastic  and  plastic  regions 
k  -  E/E' 

The  tangential  strain  at  the  elasto-plastic  boundary  (r=R)  must  have  the  same 
value  in  both  the  elastic  and  plastic  regions.  Since  the  values  of  the  Young's 
modulus  are  different  in  the  two  regions,  there  will  be  a  discontinuity  in  the 
temgential  stress  at  r  =  R. 

Let  be  the  radial  stress  at  the  elasto-plastic  boundary  (r=R) 

The  circumferential  strain  at  r  =  R  due  to  stresses  in  the  elastic 

zone  is 
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The  corresponding  value  at  r  =  R  due  to  stresses  in  the  elastic  region 


is 


'GR 


i  Po  -  -  ^'^R 

ii.  ^ 


Equating  these  two  expressions  for  e 


eR 


2  p  -  k  a 
^c  u 

~  k  fN,  -  u')  +  (1  +  u) 

V 


(lt6c) 


According  to  Eq.  (UO)  Op  is  also  given  by 


=  ^Pi  (I) 


N,-l 


-  T 


(hOc) 


Equating  Eqs.  (l46c)  and  (UOc)  and  solving  for  R/a; 


R 

a 


^  Pq  - 
k  CN^  -  u')  +  (1  +  y) 


+  T) 


1/N  -1 


Ci*8c) 


Eq.  (U8o)  thus  gives  the  radius  of  the  elasto-plastic  boundary. 


For  r  >  R  (elastic  zone) 


^r  =  Po  -  ^Po  - 


^Rl^ 


(19) 


=  Po  (Po  " 


"r  =  E  ^°r  "  ^‘^0^ 


(20) 

(6) 


^0  =  E  ^°0  ■ 


(5) 
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For  a  ^  r  ^  R  (plastic  zone) 


a 

r 


(Pi  +  T) 


(1*0) 


a 


e 


N, 


a  +0 
r  u 


(33) 


The  circumferential  strain  e.  in  the  plastic  region  can  he  obtained  from  S  which 
is  given  by  Eq.  (52): 


+  T 


ir+1 

9 


(52) 


e 


0 


J, 

E' 


(53c) 


Graphical  Solution 

The  solution  to  Case  U  cem  also  be  obtained  using  the  chart.  In  this  case, 
since  the  elastic  and  plastic  zones  are  assumed  to  have  different  material  con¬ 
stants,  they  will  be  treated  as  two  different  media  arranged  in  a  concentric 
manner.  For  each  medium  (in  this  case  for  each  zone)  the  appropriate  chart  will 
have  to  be  used. 

Using  the  given  data,  the  initial  starting  point  is  located  in  the  chart  of 
one  medium.  Then  the  radial  stress  and  the  tangential  strain  at  the  boundauy 
between  the  two  media  are  determined  which  are  then  used  to  locate  the  initial 
starting  point  in  the  second  medium.  The  same  procedure  can  be  repeated  if  there 
are  more  than  two  different  media  to  deal  with. 

If  the  circuraferertial  strain  and  the  radial  stress  are  known  at  any 


radius  r  in  the  medium,  the  procedure  described  above  can  be  used  directly  to 


estimate  the  free  field  pressure  or  the  modulus  of  the  medium  required  to  limit 
the  circixmferential  strain  at  the  tunnel  wall  to  a  tolerable  value  vmder  a  given 
value  of  p^.  But  very  often  the  data  may  be  insufficient  to  locate  the  Initial 
point  in  the  chart  for  a  graphical  analysis,  as  for  example,  when  it  is  required 
to  estimate  the  circumferential  strain  at  r  =  a  in  Case  1*  for  a  given  value  of  p^. 
Under  these  conditions,  it  becomes  necessary  to  use  a  successive  approximation 
procedure.  For  the  example  cited,  sin  initial  value  of  circumferential  strain  at 
r  =  a  is  eissumed  and  the  graphical  analysis  is  carried  outward  from  the  innermost 
radius  a  and  the  vsilue  of  p^  is  estimated.  If  this  vadue  agrees  with  the  value 
given  in  the  problem,  the  estimate  of  the  circumferential  strain  at  r  =  a  is 
correct.  If  not,  a  revised  value  of  the  circimiferential  strain  is  assumed  and  the 
procedure  is  repeated  until  the  assumed  circumferential  strain  is  consistent  with 
the  given  data. 

Example  5 

Data: 


Unlined  Tunnel:  a 

= 

8’ 

Po 

= 

2000  psi 

E 

= 

6x10^  ps: 

E' 

= 

3x10^  ps: 

y 

= 

1/3 

y' 

= 

1/3 

= 

1200  psi 

_ 

U 

<P 

Pi 

= 

0 

30 


Required!  =7 


Solution: 


T  =  =  UOQ  psi 


=  E/E'  »  2 


R  (  1*000  -  21*00  .  /,.„ 

I  ■  \[2  X  3-/3- 1  1/3  ^  '•“J /  '•“J 


(l+8c) 


=  1.136 


(Compare  with  the  value  1.3l*  obtained  for  Example  2  with  E  =  E'  and  constant), 


At  r  =  u 


S  3 


Tioo  =  8 


3  +  5  (1 


S  =  251*0  psi. 


^  [251*0  +  2/3  X  Cj  =  8.5  X  10"** 

®  3x10^ 


(Compare  with  the  value  IbxlO*  obtained  for  Example  2  with  E  =  E'  and  cr^  = 
constant ) . 


31 


32 


‘eR-  <1-*'''  “r 


or 


'eR 


"  F"  Sr  +  (1  -  m')  Or 


3x10 


1755  +  ^  X  185 


=  626  X  10 


-6 


Now  in  the  elastic  region 


S  ^  "eR  -  -^R 


=  6x10^  X  626x10“^  -  I  X  185 


=  (3756  -  123)  10  =  3633  ’  ai 


2p^  =  Sj^  +  Op  =  4003  psi 


=  2001.5  =  2000  psi  (given) 


This  means  that  the  initial  assumed  value  of  e 


is  correct. 


r=a 


If  the  value  of  p  obtained  for  an  assumed  value  of  Go 
O  G 


does  not 


r=a 


agree  with  the  given  value  of  p^,  the  procedure  has  to  he  repeated  with  a  second 


trial  value  for  e 


.  This  repetitive  procedure  is  continued  until  aill  the 


I  r=a 

input  data  are  satisfied. 
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Case  5.  In  this  case  the  tunnel  Is  assiuued  to  be  initially  surrounded  vith  a 
circuleg*  zone  of  loosened  material  heiving  elastic  constants  which  are 


different  from  those  of  the  Insltu  material. 

The  case  could  represent  a  tunnel  with  a  "destressed  zone"  from  static 
loading  which  is  subsequently  loaded  with  another  free-field  stress  such  as  a 
dynamic  loading  over  a  large  area. 


Let 


a  =  radius  of  tunnel  opening 

R  =  radius  of  elasto-plastic  boundary 

=  radius  of  loosened  zone 

=  radius  of  elasto-plastic  boundary  in  the  loosened  medium,  if 

any,  when  R  > 

=  uniform  free-field  pressure 

p^.  =  interned  pressure  at  r  =  a 

E  =  Young's  modulus  of  intact  medim 

E'  =  YoTing's  modulus  of  loosened  medim 

M  =  Poisson's  ratio  of  intact  medium. 

y'  =  Poisson's  ratio  of  loosened  medium, 

=  unconfined  compressive  strength  of  the  medium  (assumed  to  be 

the  same  in  both  the  intact  and  loosened  zones) 

Ij  =  1  +  sin  (p 

(p  1  -  sin  ()) 

4)  =  angle  of  shearing  resistance  of  the  medixu^'  (assumed  to  be 

the  same  in  both  the  intact  and  loosened  zones) 


In  the  analysis  of  t}..e  ;  r'esent  problem,  two  distinct  cases  have  to  be 
considered: 

a)  R  R^  b)  R  >  R^ 
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Case  (a)  R  i  Rl 


The  elastic  zone  (r  ^  R)  in  this  case  consists  of  two  rings  of  material 
with  two  different  values  of  Young's  modulus.  Therefore  the  stresses  and  strains 
in  the  elastic  zone  cam  be  computed  using  the  elastic  solution  obtained  by 
Savin  (l96l)‘. 

For  R  <  r  <  Ri 

_  —  — »  A 


^r  =  ^  (Po  -  “r)  ^1 


(61) 

(62) 


where 


C.  = 


[C1  .  x^)  (n"  -  RiVr=)]/l> 


=  [(1  +  Xg) 


(63) 


iSU) 


and 


D  = 

2  (a  -  1)  -  n^  jja  -  l)  -  (l  +  a)j 

(65) 

3^  Kl 

1  +  m' 

(66) 

Xg  = 

1  +  p 

(67) 

n  = 

R^/R 

(68) 

a  = 

E  (1  +  u') 

E'  (l  +T) 

(69) 

°R  = 

"r| 

r=R 
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For  r  ^  R-j^ 


=  Po  - 

=  Po  +  CPq  -  cJj^)  C3  (71) 


where 

C  S= 

(n 

1 

^  -  1)  (1  +  x^) 

(72) 

D 

2 

r 

The  strains  in  the  elastic  zone  can  be  obtained  from  the 

stresses  as 

follows 

For  R  _< 

r  ^  P.^ 

G  = 

r 

1  t  t 

F  ^°r  ■  P 

'e^ 

(73) 

11 

<x> 

(Oq  -  P'  0 

'r^ 

(71*) 

For  r  _> 

e  = 

r 

1  ‘“r  -  "  '=e> 

(75) 

"e  = 

E  -  P  “r^ 

(76) 

However,  in  Eqs. 

(61), 

(62),  (70)  and  (71)*  Op  is  still  an  urJcnown  sind 

therefore  has  to  be  obtained  before  the  stresses  and  strains  can  be  determined  in 
the  elastic  region.  This  can  be  done  by  equating  the  circumferential  strains  at 
r  =  R  in  both  the  elastic  and  plastic  regions.  Thus 


‘“eB.  -  “r>  = 


__1 

E' 


JL 

E' 


(77) 


where 


la  the  elastic  region 


“^eRe  " 


r=R 


°eRp  ■ 


r=R 


In  the  plastic  region 


Eq.  (77)  gives 


By  Eq.  (62)  +  (p^  -  Oj,) 


r=R 


=  °R  *  ^Po  -  ^ 


(78) 


where 


C  =  (1  +  X2)  •  2  n  /D 


In  the  plastic  region 


"ORp  ■  "r  •  % 


CT91 


Equating  Eqs ,  (78)  and  (79)  1  aJid  simplifying 


C  p  -  o 

u 


N.  -  1  +  C 


(80) 


Eq.  (80)  can  now  be  used  in  conjunction  with  Eqs.  (61)  throxigh  (76)  to  define 
completely  the  stresses  and  strains  in  the  elastic  zone. 

From  Eq.  (Uo)  it  is  known  that 


=  (Pp  ♦  T)  [l] 


N.-l 


(81) 


Eqs.  (80)  and  (81)  can  now  be  combined  to  yield 
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The  value  of  R/a  obtained  from  Eq.  (82)  should  be  less  than  R^/a  for  the  above 
analysis  to  be  valid.  The  analysis  of  the  problem  when  R  >  is  considered 
under  Case  (b) . 

The  stresses  in  the  plastic  zone  (a  ^  r  ^  R)  are  given  by  Eqs.  (Uo)  and 
(33).  The  circumferential  strains  in  the  plastic  zone  can  be  determined  from 
Eq.  (52)  and 


S  =  E'  Cg  -  (1  -  y)  0^ 


(83) 


Case  (b)  R  >  R]^ 

When  R  >  R^,  the  loosened  zone  may  exist  in  any  one  of  the  following 
three  states,  depending  on  the  magnitude  of  the  radial  stress  at  r  =  R^: 

(i)  Completely  plastic 

(ii)  Partly  plastic  and  peirtly  elastic 

(iii)  Completely  elastic 

Therefore  it  is  necessary  to  first  determine  the  magnio’’.de  of  The  procedure 

for  determining  o^^  is  given  below. 

Rased  on  Eq.  (48)  it  may  be  written  that 


__2 _ 

N.  +  1 


Pq 

"ri 


(84) 


Similarly  from  Eq.  (52)  it  can  be  derived  that 


38 


'HI 


2N. 


2N, 


(N  -  1)  +  CN.  +  1) 

9  9 


(85) 


vhere  -  E  <Jh1 


(86) 


€Uld  C  =  £ 

^eRi 


r«Rl 


Eqs.  (81t),  (85)  and  (86)  oeui  be  combined  to  give 


^eRl 


1 

E 


2N^' 

Na-1 

N  -1 

r  2  Po-^T  • 

y 

*  2N 

9 

(K  -1)  +  (N  +1) 

<P  9 

N  tl  • 

- 

. 

K 

I 


+  (1-w)  o 


R1 


(87) 


As  a  genereu  case,  let  it  be  assumed  that  the  loosened  zone  becomes  peirtly 
plastic  and  partly  elastic.  If  R^  is  xhe  radius  of  the  elasto-plastic  boundary 
in  the  loosened  zone  the  following  relationships  must  be  satisfied. 


N  -1 


°R2  '  '^r 


=  (p,  +  T) 


r=R, 


R. 


-  T 


(UOd) 
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r 


in  the  elastic  loosened  zone 


0 

eRle 


'^eRl  ~  E'  ^°0Rle  “  °R1^ 

Equating  Eqs.  (8?)  and  (5d) 

‘^eRle  "  °R1  ■  ~  ^ 

It  can  also  be  shown  that 


r=R, 


2  2 

Rl  +  Rg 

R  2  2 

^1  -  ^2 


-  a. 


2R 


R  2  R  ^ 
"l  "  ^2 


(88) 


{5d) 


(89) 


(N. 


'’E2  *  “a 


=  2 


^R1  *  ~  °P2 

Hl^  -  h/ 


(90) 


Eqs.  (89),  (90)  and  ( l*0d)  can  now  be  solved  to  yield  the  values  of  the  two 
unknowns  and  R^. 

If  the  loosened  zone  remains  entirely  elastic  or  entirely  plastic, 
Eqs.  (l*0d),  (88)  and  (90)  are  not  veilid  and  cannot  be  used  to  determine  Om  • 

Ill 

When  the  loosened  zone  is  entirely  el8u?tic 


0 


0Rle 


.  a^)  - 


(91) 


This  value  of  cein  now  be  substituted  into  Eq.  (89)  to  get  the  value  of 

Substitution  of  in  Eq.  (8^4)  gives  the  value  of  R. 
ill 


When  the  loosened  zone  is  entirely  plastic,  the  values  of  R,  o_  and 
cam  be  obtained  directly  from  the  following  equations; 


2  Po  - 


»  -H  p,  »  T 


(2  /  (H^  .  1) 


The  value  of  R  obtained  from  the  preceeiing  equations  should  be  greater 
them  R^.  If  not,  the  procedure  given  vmder  case  (a)  has  to  be  used. 


Stresses  amd  Strains 


r  >  R 


Or  =  Pq  -  (Pq  -  Op)  - 


«e  =  Po  ^Po  -  '^R^  'r 


S  =  E 


where 


Or  -  0^ 


(2  p  -  o  ) 
•^o  u 

N,  +  1 


Ul 


I  ^ 


<  r  <  R 


2  (P^  +  T)  , 

°r  “  N,  +  1  R  " 


a,  =  o  •  N,  +  o 
G  r  <j)  u 


^  +  T  I  R 


(»♦  *  i: 


I 

+  (N.  +  1) 

<P 


wh  .re 


S  =  E  Cq  -  (1  -  y) 


a  <  r  c  R-| 


0  =  (p.  +  T)  -  -  T 

r  [&J 


(100) 


0.  =  a  •  N,  +  a 

e  r  if  u 


(N  -1)2  (p.  4.  T)  r  2H 


^"ri  ^1 


(101) 


where 


E'  Eg  -  (1  -  y') 


®R1  ~  ‘^GRl  ‘^Rl 


and 


^eRl  ”  ^6' 


aa  obtained  from  Eq.  (99). 


r=Rl 


Partly  plastic  and  partly  elastic 


^2  ^  i  % 


°R1  •  ^1^  -  °R2  •  ^2  ^‘^Rl  -  ^1^  ^2 


2  2 

Rl  -r/ 


(102) 


‘^Rl  •  ^1^  -  ^R2  •  ^2  ,  ^°R1  -  ^RP^  ^1^ 

*^0  ~  „  2  „  2  2  ^  2s  2 


\  -  «2 


(R^"  -  R^^)  r' 


(103) 


"r  *  F  <°r  - 


(lOlf) 


=  F  (cTq  -  P'  O^) 


(105) 


a  ^  r  j5_  R2 


(p  +  T)  f|  -  T 
1  (_aj 


(100) 


0-,  =  0  •  N,  +  a 

6  r  <()  u 


(98) 


N,+l 


Pi>T 


_i _ 

2N. 


(».  -  1) 


2N. 

w 

laJ 


2N, 


(N,  +  1) 


(106) 


S  --  E’  Cg  -  (1  -  p') 


Us 


where 


Fully  elastic 


a  = 

r 

A  -  B/r^ 

(107) 

A  +  B/r^ 

(108) 

where 

o  2  2 

A  — 

''ri  *  ^1  -  Pi  *  ^ 

(109) 

A,  — 

D  2  2 

^1  -  ^ 

(Opi  -  p. )  R^2  ^2 

(no) 

o  — 

(r/  -  a^) 

c  = 

r 

|t 

(10b) 

F-  <»0  -  "r* 

(105) 

The  preceeding  peiragraphs  thuii  give  a 

complete  and  closed  analytical 

solution  for  Case  5* 

This  problem  can,  however,  be  solved  more  readily  by  tne 

graphical  procedirre  using  appropriate  charts. 

The  details  of  the  procedure  are 

the  same  as  explained  at  the  end  of  Case  U. 

Example 

The  input  parameters  for  this  example  are  the  same  as  for  Example  1  except 
that  the  tunnel  opening  is  surrounded  by  a  loosened  zone  8  feet  thick  having 

bit 


the  following  properties : 


V. 


E'  =  2x10^  psi 

U'  =  1/3 

0  =  2000  pai 

u  ^ 

tr  =  h 

9 

=  ? 

r=a 

Solution: 

This  problem  can  be  solved  analytically  using  the  relationships  derived 
for  Case  5.  But  the  procedure  is  long  and  time-consuming.  Therefore  the  easier 
graphical  proced\ire  is  used  in  this  example. 

As  has  been  explained  earlier,  the  solution  to  this  problem  can  only  be 
obtained  by  successive  trials,  the  number  of  trials  being  dependent  on  how  close 
the  first  trial  value  is  to  the  actual  value.  In  the  present  case  let  us  assume 

eg  =  0.05 
r=a 


Then  at  r  =  a 

% 

2000  psi 

^0 

s 

0.05 

0 

0 

r 

E 

= 

E'  =  2x10^  psi 

V  ' 

= 

U  =  1/3 

* 

«  • 

5 

= 

2x10^  X  0.05  -  0  =  10^  psi 

(9) 

P 

= 

2000  +0  =  2000  psi 

(3U) 

For  these  values  of  S 

and  P, 

read  from  Chart  (II,  =  U) 
w 

At  r 


For  Q 


At  r  = 


At  r  = 


Q  =  8600  pai 


-  =  1.725 


=  2a  (=  16') 


^  =  0.063 

^1 


8600  psi  and  — 
r 


R  R 


—  =  0.863,  read  from  the  chart 

h 


(loosened  medium)  =  7700  psi 


r=R, 


<  1  ,  the  material  is  elastic  at  r  =  R, 


r=K, 


2  -  S  2  (Q  -  T)  -  S 


=  (Q  -  T)  -  I  (19) 


(8600  -  667)  -  =  U083 


psi . 


in  the  loosened  zone. 


S  +  (1  -  li') 


M  )  0 


2x10 


^  [7700  +  I  X  lt083j 


=  5211  X  10 


-6 


Rj^  in  the  unloosened  zone 


1*6 


(9) 


S  =  E  Eg  -  (1  -  u)  0^ 

=  6x10^  X  5211  X  lO"^  -  I  X  1*083 
=  285UU  psi 

P  =  a  +  (N.  -  1)  0  =  2000  +  3  X  1*083  (3l*) 

U  9  T 

=  11*250  P3i 

For  these  values  of  S  and  P,  the  chart  for  N.  ==  U  gives 

<(> 

Q  =  16900  or  p^  =  Q  -  T  =  16900  -  667  =  16233  psi 

This  value  is  sufficiently  close  to  the  given  value  of  p^  of  161*00  psi.  So  the 

assumed  value  of  =  0.05  is  very  near  the  correct  value  and  the  calcula- 

r=a 

tions  need  not  be  repeats • 

Case  6.  Lined  tunnel  vith  a  circular  zone  of  loosened  material. 

A  tunnel  may  be  provided  with  a  liner  to  limit  the  strains  at  the  inner 
surface  of  the  tunnel  within  allowable  limits.  The  linings  may  be  of  concrete 
or  steel  or  a  combination  of  both.  In  this  analysis  the  liner  materials  also  will 
be  assumed  to  behave  elasto-plastically .  The  analytical  solution  of  the  stresses 
and  strains  in  a  tunnel  system  consisting  of  one  or  more  sets  of  liners  and  a 
loosened  zone  of  material  surroxinding  the  liners,  is  very  complex,  lengthy  and 
cumbersome.  However,  the  graphical  solution,  using  the  charts  of  the  type 
presented  earlier,  offers  a  relatively  simple  means  of  solving  the  problem. 

J*7 


The  analysis  must  proceed  outward  from  an  inner  element  or  inward  from  an  outer 


element  of  the  tunnel  system.  For  each  element  of  the  system  the  chsurt  with  the 
appropriate  value  of  F,  has  to  be  used.  At  least  two  quantities! such  as  the 
circumferential  strain  or  the  radial  stress,  are  required  to  locate  the  initial 
starting  point  for  the  analysis  of  any  one  element.  Once  these  quantities  are 
known,  they  can  be  used  to  estimate  the  values  of  the  radial  stress  and  circum¬ 
ferential  strain  at  the  boundary  with  the  next  element,  which  are  in  turn  used 
to  locate  the  starting  point  for  the  analysis  of  that  next  element.  This  process 
can  be  repeated  until  the  stresses  and  strains  in  the  whole  tunnel  system  are 
known.  The  following  example  is  worked  out  to  illustrate  the  approach  mentioned 
above . 

Example  ^ 

The  dimensions  and  rock  properties  assumed  in  this  example  (Fig.  8)  are 
the  same  as  for  Example  U  except  that  the  tunnel  opening  is  provided  with  a  12" 
thick  liner  of  concrete  having  the  following  properties: 


E  = 
c 

4x10^  psi 

u  = 
c 

1/3 

0  = 
u 

'yOOO  psi 

IJ.  = 

1* 

Eeq  aired: 

e  !  -9 

0 

1  r=a 

U8 


Solution; 


Let  us  consider  the  concrete  liner  first.  Let  us  assuue,  as  a  first 
trial  value, 


G 


e 


0.03 


It  is  known  that  at  r  =  7' 


S  =  U  X  10^  X  0.03  -  0  =  1.2  X  10^  psi  C9) 

P  =  5000  +  0  =  5000  psi  (3k) 

For  these  vedues  of  S  and  P,  read  from  chart  IN, (concrete)  = 

9 

Q  =  16500  psi 

Y  =  1.57,  thus  R  =  11' 

At  r  =  8'  £  =  =  1.375 

For  Q  =  16500  psi  and  ^  =  1,375»  read 

S  =  6^000  psi 

P  =  7600  psi  =  5000  +  3  (3i*) 

Therefore  at  r  =  8' 

.  .  =  867  psi 


k9 


[6U000  +  2/3  X  867] 


=  16114?  X  10" 


Let  us  now  consider  the  loosenei  zone. 


At  r  =  8' 


Cq  =  I61U7  X  10" 


=  867  psi 


.  S  =  2  X  10^  X  l6llt7  X  10“^  -  2/3  x  86? 


=  31720  psi 


P  =  2000  +  3  X  867  =  1*600  psi 


For  these  values  of  S  and  P,  read  from  the  N  =  1*  chart  for  the  loosened  zone, 


9250  psi 


=  I.3U 


R  ^  10.72' 


At  r  =  16' 


10.72 


=  0.67  <  1 


For  Q  = 


9250  psi  and  p 


=  0.67,  read 


S  =  5000  psi 


As  —  ”  C.67  <  1,  the  value  of  P  is  not  valid,  but  can  be  calculated  using 

tlie  relation 


a  =  P  -  cr  wnere  p 
r  -^o  2  o 


=  Q  -  T 


50 


r 


Thus  at  r  =  l6* 


-  T)  -  I  *  (9250  -  667)  -  -^222. 


6083  psi 


*^0 


1  r 


|5000 


+  4  X  6083 


=  lt528  X  10* 


Let  us  now  consider  the  intact  (non  loosened)  medltun. 


At  r  =  16' 


Cq  »  i*528  X  10* 


=  6083  psi 


.  .  S  =  6  X  10^  X  1*528  X  10"^  "  3  ^ 


=  23113  psi 


P  =  2000  +  3  X  6083 


=  20250  psi 


For  these  values  of  S  and  P, 


Q  =  18200  psi 


=  Q  _  T  =  18200  -  667  =  17533  psi 


i  j 


But  the  veLLue  of  p 

assumed  value  of  £„ 
y 


specified  in  the  problem  is  161400  psi.  This  means  that  th 
=  0.03  is  slightly  on  the  high  side. 

r=7' 


So  let  us  now  assume  that  at  r  =  7' 

Eg  =  0.028 


A  similar  aneilysis  gives  the  following  values: 

At  r  =  7*  In  concrete 

S  =  1.12  X  10^  psi  P  =  5000  psi 

At  r  =  8'  In  concrete 

4 

S  =  5.95  X  10  P  =  7500  psi 

Op  =  833  psi 

Cg  =  150114  X  10"^ 


In  loosened  zone 

S  =  29I472  psi  P  =  I45OO  psi 

At  r  =  16'  In  loosened  zone 

S  =  l4.7xl0^psi 
0^  =  5883  psi 

e-  =  I43II  X  10"^ 

y 


52 


In  intact  zone 


S  ■  psi  P  =  19650  pai 

=  16433  (“  16400)  pal 

Thus  the  assumed  value  of  Cg  (at  r  =  7')  “  0.028  is  very  neao*  the  actual 

value. 


Note;  For  relatively  thin  steel  liners,  the  analysis  is  slightly  simplified 
because  the  radial  pressure  exerted  by  the  thin  steel  liner  against  the  mediim 
next  to  it  is  given  by 


ra 


=  —  E 
a  a 


or 


a 

y 


whichever  is  less.  In  the  above  equation 


a  =»  radial  pressure  exerted  by  steel  against  the 
adjacent  mediuni 

h  =  thickness  of  steel  lining 
8 

Oy  =  yield  stress  of  steel  lining 
a  =  radius  to  outside  of  steel  lining 
Eg  =  Young's  modulus  of  steel  lining 
e.  =  circujnferential  strain  of  the  steel  liner 

U 


Case  T.  Tui.inel  provided  with  back  packing:  Loosened  zone  present. 

This  case  is  very  similar  to  Case  5,  the  only  difference  being  that  the 
radial  pressiure  p^  at  r  =  a  is  #  0.  The  back  packing  material  usually  has  a 
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very  low  yield  value  and  serves  to  exert  an  eq.ual  all  around  pressure  on  the 
inner  surface  of  the  tunnel  over  a  considerable  range  of  strain.  The  analysis 
for  the  present  case  is  illustrated  by  the  following  example  problem  and  is  very 
similar  to  that  of  Example  U. 


Example  6 


The  input  parameters  are  the  same  as  those 
radial  pressure  due  to  back  packing  is  I50  psi  (p^ 


of  Example  b  except  that  the 
=  150  psi) . 


Required: 


e 


6 


r-a 


<7 


Solution: 

Assume  at  r  =  a  =  8' ,  ^6  " 

Then  at  r  =  8' 

Z  ^  ?  X  10  X  .04  -  ^  X  150  =  79900  psi  (9) 


P  =  2000  +  3  X  150  =  2450  psi  (34) 

For  these  values  of  S  and  P  read  from  chart  (tl  =4) 

•t 

Q  =  9000  psi  =  1.64* 

R  =  13.12' 


At  r  =  16' 


R  ^  13.12 

r  1^.0 


0.82  <  1 
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For  Q  ■  9000  psi,  ^  -  0.82 


S  -  7250  pal 


-  (Q  -  T)  -  I 


(9000  -667)  -  3625 


1*708  pal 


2x10 


b  *  I  -  H 


5195  X  10 


-6 


Now  for  the  intact  material  ct  r  »  I6' 


*  1*708  pal 


Cg  »  5195  X  10 


-6 


S  =  6  X  10^  X  5195  >•  10"^  "  3  * 


(9) 


=  28040  pai 


P  =  2000  +  3  X  1*708  »  16121*  pal 


C3I*) 


For  these  veLlues  of  S  and  P 


Q  »  18000  pai 


a  Q  -  T  =  17333  psi  Cin  compariaon  with  the  actual  value  of 

l61*00  pai) 


55 


Therefore  the  assumed  value  of  at  r  =  a  =  8’  has  to  "be  revised.  So  let  us  now 
assume  at  r  =  8*  ^0  “  0.038. 

A  similar  analysis  leads  to  a  value  of  =  16633  psi.  So  by  a  slight 
extrapolation  it  can  be  found  that  for  p^  =  l6U00  psi, 

Eq  «  0.0373  or  3.7/^. 

r=8' 
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I  iju 


Chapter  3 
Concluaiona 


General 


The  analysis  presented  in  the  previous  chapter  yields  a  means  whereby 
stresses  and  strains  can  be  determined  £a*oijnd  a  circular  txmnel  in  a  Coulomb* 
Navler  material  idiich  increases  in  volume  at  failure.  This  analysis  may  be  losed 
to  design  liners  if  the  insitu  rock  properties  are  known.  The  results  of  the 
analysis  are  dependent  on  the  rock  properties  assumed;  therefore,  recommenda- 
tions  are  given  in  this  chapter  on  the  selection  of  the  approrlate  rock  mass 
properties  to  use  in  this  analysis.  The  effects  of  a  non-uniform  system  of 
free-field  stresses  are  also  discussed. 

Shear  Strength  Properties 

The  values  of  o  and  N.  to  be  \ised  in  the  analysis  ore  related  to  the 
u  9 

properties  of  the  Jointed  rock  mass  sxurrounding  the  opening.  is  given  by 

1  "  ~si^^  where  the  appropriate  angle  of  frictional  resistance  should  be  taken 

as  the  angle  of  frictional  resistemce  along  the  Joints  or  discontinuities  not 
the  angle  of  internal  friction  derived  from  triaxiol  tests  on  intact  samples  of 
rock.  The  value  of  the  unconx-ined  compressive  strength  of  the  rock  mass,  is 
a  function  of  the  ratio  of  tunnel  diameter  to  Joint  spacing  as  shown  in  Fig.  9* 
The  larger  the  ratio  of  tunnel  diameter  to  Joint  spacing,,  D/S,  the  smaller  the 
value  of  appropriate  for  design.  If  D/S  is  very  small  the  value  of  can 
approach  the  unconfined  strength,  q^,  of  intact  entitles  of  the  rock  surrounding 
the  tunnel.  Similarly  if  D/S  is  very  large  approaches  zero  and  the  ahear 
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strength  of  the  rock  mass  approaches  the  sheeir  strength  along  the  Joints.  Fig. 
10  gives  a  straight  line  relationship  hetween  the  ratio  of  and  the  ratio 

of  d/s  which  can  be  used  to  select  a  value  of  for  design  if  and  D/S  are 
known.  Tliis  relationship  is  based  on  field  experience  in  locations  where  dis- 
placenent  and  strain  measurements  have  been  made. 

Elastic  Properties 

The  method  of  analysis  presented  in  Chapter  2  is  also  quite  dependent  on 
the  "effective"  modulus  of  elasticity,  E,  selected  for  the  rock  mass.  A  method 
for  selecting  the  deformation  modulus  of  a  rock  mass  has  been  given  by  Deere, 
Hendron,  Patton  and  Cording  (1967).  By  this  method  the  rock  quality  of  the  rock 
mass  must  first  be  assessed  quantitatively  in  terms  of  the  Rock  QueLLity  Designa¬ 
tion  (RQD)  or  the  Velocity  Ratio  as  described  by  Deere  et  al  (1967).  After  the 
rock  quality  has  been  determined  Fig.  11  should  be  entered  on  the  abscissa  using 
the  RQD  value  or  the  square  of  the  velocity  ratio  and  the  reduction  factor 
E  /E  .  should  be  determined  from  the  dotted  line  shown  in  Fig.  11.  The 

X*  S6XS 

reduction  factor  is  the  ratio  of  the  deformation  modulus  of  the  rock  mass  E 

r 

to  the  dynamic  value  of  Young's  modxilus  calculated  from  P  wave  velocities 

measured  in  seismic  surveys.  The  deformation  modulus  of  the  Jointed  rock  mass 

can  then  be  taken  as  the  product  of  the  reduction  factor  and  E  .  . 

^  seis 

This  analysis  is  not  very  sensitive  to  the  value  of  Poisson's  ratio 

selected,  which  is  fortunate  because  there  is  little  known  about  the  selection 

/ 

of  a  Poisson's  ratio  for  an  insitu  rock  mass.  It  is  recommended  that  a  Poisson's 

:!ati '  of  o'.'  ut  0.3  be  jf.ed. 
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Dllatancy 


DiJ.atancy  of  the  rock  mass  at  failure  is  accounted  for  in  the  elasto- 
plastic  analysis  by  the  normality  condition  expressed  by  Eq.  2h.  E^erlmental 
evidence  hovever  siiggests  that  the  increases  in  volume  resulting  from  the 
normality  condition  are  too  large  compared  to  the  behavior  of  real  rocks.  Thus 
the  radial  displacements  or  tangential,  strains  computed  by  this  theory  acre 
a  conservative  upper  boimd.  On  the  other  hand  the  chaxts  developed  by  Newmairk 
(1969)  sore  based  upon  no  volume  chsuxge  due  to  the  plastic  strains  at  failure 
and  give  radial  displacements  and  tauigentlal  strains  which  are  too  small  and. 
should  be  taken  as  a  lower  bound.  The  cheu*t3  given  in  Figs.  3-7  would  be 
identical  to  Newmark's  chsirts  if  the  lines  of  constant  Q  in  the  elastic  region 
axe  extended  as  straight  lines  into  the  plastic  region.  The  proper  aaoxint  of 
dllatancy  to  include  in  calculating  strains  and  displacements  is  somewhere 
between  that  given  by  Newmark  (1969)  and  the  anedysis  given  herein. 

Non-uniform  Stress  Conditions 

For  real  design  problems  the  free-field  principal  stresses  in  a  plane 
perpendlculsu:  to  the  tunnel  axis  may  be  o^  euid  where  is  the  ms^or  principal 
free-field  stress  and  is  the  minor  free-field  principal  stress.  The  solution 
for  this  problem  cannot  be  obtained  in  closed  form  as  was  done  in  this  paper 
for  the  case  of  ^3^  “  “^3  ~  Pq*  Reyes  (1966)  has  solved  several  problems  for 
unllned  openings  in  a  medium  with  Coulomb-Navi er  failure  properties  and 
dllatancy,  as  described  in  this  report,  for  non-uniform  free-field  stress  con¬ 
ditions.  A  conparison  of  the  uniform,  free-field  stress  solution  presented 
herein  with  the  solutions  presented  by  Reyes  has  shown  that  the  maximum  dlametred 
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strain  6^/r  (max)  across  the  tunnel  in  the  non-uniform  stress  field  (free- field 

stresses  =  and  o^)  is  closely  approximated  by  the  solution  given,  herein  if 

p^  is  assumed  to  be  equal  to  the  me^Jor  principal  free- field  stress,  The 

distortion  of  the  tunnel  from  a  circular  shape  can  also  be  approximated  within 

20  %  if  the  minimum  diametral  strain  6  /r  (min)  is  teiken  as 

r 

/r  ^  /r 

r  (min)  r  (max) 

Since  the  ratio  of  protective  structures  problems  may  range  from 

about  1/3  to  2/3  the  minimum  diametral  strain  of  the  tunnel  may  be  approximated 
for  preliminary  design  purposes  as  being  about  1/3  to  2/3  the  maximum  diametral 
strain  across  the  tunnel.  Reyes  solution  has  also  shown  that  the  maximum 
diametral  strain  of  the  tunnel  occiars  across  a  diameter  parallel  to  the 
direction  of  and  the  minimum  diametral  strain  occurs  across  a  diameter 
parallel  to  Oy 

For  preliminary  design  piarposes  it  is  felt  that  the  approximations 
given  in  this  section  for  estimating  tunnel  distortions  for  the  non-uniform 
stress  field  are  at  least  as  accurate  els  the  initial  assumption  of  the  ratio 
of  ^appropriate  for  design  problems  in  protective  structures.  If  more 

accurate  estimates  are  desired,  then  time  consuming  and  expensive  finite 
element  calculations  similar  to  those  used  by  Reyes  (1965)  must  be  performed. 
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FIG.  I  CIRCULAR  TUNNEL  SUBJECTED  TO  AN  ISOTROPIC  ERFF-FiFl  O  RTRT 


FIG.  2  ELASTIC- PLASTIC  STRESS  DISTRIBUTION  AROUND  A  CIRCULAR  TUNNEL  IN  SAND 


FIG.  4  GRAPHICAL  ILLUSTRATION  OF  RELATIONS  AMONG  P,  S,  AND  Q  FOR  N^=3 
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FIG,  5  GRAPHICAL  ILLUSTRATION  OF  RELATIONS  AMONG  P,  S,  AND  0  FOR  N^=4 


. ‘ ‘WeilBMfr 


iK4 


_ 

..  .jm  .UlIttMIfUUUUSR’^L 


Si 

JPj 

Ml4kM| 


! 


ililr 


mt 


lilM’.ltNKl 


mill! 


mi 


M\ 


& 


M 

i 


*  ^  6  7  8  9  KD.OOO 


4  5^789  CO.OOC' 


i  4  *1  8  ’  ^JOO,OO0 


S-Ztg  -  (  I 
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FIG.  8  GEOMETRY  OF  PROBLEM  ASSUMED  IN  EXAMPLE  5 


Cl 

•e- 


5  c 

*1  8 
or  c 
LU  r) 


<u  ^ 


C 

8 

0>  C 


♦-  O 

c  q: 


c  c 


o  o  o 
li:  K  -I 


I—  LJt3 

LU  |sj  3 

LJS  ii. 

S  c/)< 

^  v^a 

Qos 

_l  Oi 


HI 

ii 

m|co 


2  <-> 

“  ^ 
CO  CL 

2  CO 


o  z 

Q  ^ 

^  O 
q:  h- 


71 


Reduction  Factor 


10 


■■■ 


Rock  Quality,  (Vp/V|_),  ROD 
DYNAMIC  TEST  DATA 


CEDAR  CITY  TONALITE 

QUARTZ  MONZONITE  -  CLIMAX  STOCK 

BUCK BOARD  MESA  BASALT 


FIG.  M  RELATION  BETWEEN  ROCK  QUALITY  AND  REDUCTION 
FACTOR 


